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Modeling  and  Characterization  of  Phonon  Transmission  and  Generation 
across  Engineered  and  Strained  Interfaces  for  Developing  Structure-Property 

Relations  of  Functional  Nanostructures 

(FA9550-1 1-1-0109) 

Performance  Period:  June  15,  2011  -  December  14,  2014 
Ronggui  Yang 

Department  of  Mechanical  Engineering,  University  of  Colorado  at  Boulder 

The  Air  Force  is  increasingly  relying  on  the  benefits  of  novel  multifunctional  materials.  These 
materials  are  being  pushed  to  their  thermal  limits  in  progressively  more  demanding  roles. 
Therefore,  there  is  a  strong  need  to  develop  modeling  and  characterization  tools  for  better 
understanding  of  thermal  transport  in  nanostructured  functional  materials.  The  objective  of  this 
project  is  to  develop  an  integrated  modeling  and  characterization  tool  to  study  phonon 
transmission  and  generation  across  engineered  and  strained  interfaces.  Such  a  tool  will  be  critical 
in  understanding  the  fundamentals  of  energy  transport  and  in  developing  structure-property 
relations  of  functional  nanostructures  that  could  enable  novel  design  of  multi-functional 
nanostructures  with  superior  thermal  properties,  such  as  super-thermal  insulators,  high  efficiency 
thermoelectric  materials,  and  super-thermal  conductors. 

The  Innovation  and  Accomplishments  of  this  project  include: 

1) .  Developed  modeling  and  simulation  tools  for  thermal  and  thermoelectric  transport  processes 
across  interfaces  and  in  nanostructured  materials.  In  particular,  we  have  developed  a  novel 
approach  to  predict  phonon  transport  properties  across  dissimilar  interfaces.  Our  approach  is 
based  on  non-equilibrium  Green's  function  theory  (NEGF)  which  provides  a  wavelength- 
dependent  phonon  transmission  function.  This  approach  has  been  integrated  with  the  molecular 
dynamics  simulations  and  the  first-principles  simulations  to  predict  phonon  transmission 
function  across  strained  and  engineered  materials  interfaces.  The  obtained  wavelength-dependent 
phonon  transmission  function  will  be  integrated  into  the  Boltzmann  transport  model  to  study  and 
design  large-scale  multidimensional  systems  with  embedded  nanostructures.  The  results  were 
published  in  Journal  of  Physics:  Condensed  Matter  (Vol.  24,  Art  #  155302,  2012),  Physical 
Review  B  (Vol.  86,  Art  #054305,  2012),  and  a  recently  accepted  paper  in  Physical  Review  B 
(Xiaokun  Gu,  et  al.  Phonon  transmission  across  Mg2Si/Mg2Sii-xSnx  interface:  a  first-principles- 
based  atomistic  Green’s  function  study,  http://arxiv.org/abs/1501.04084.  May  2015) 

2) .  In  an  earlier  work  published  in  Nature  Materials  (Vol.  9,  pp.  26-30,  2010),  we  developed  a 
pump-and-probe  characterization  tool  using  table-top  femtosecond  extreme  ultraviolet  (EUV) 
laser  to  study  nanoscale  thermal  transport.  In  the  experiment,  an  array  of  nanowires  made  by 
electron-beam  lithography  is  heated  with  an  ultrafast  laser  pulse,  and  the  heat  flow  from  the 
nanowire  array  into  a  bulk  substrate  is  monitored  by  diffracting  ultrafast  coherent  EUV  beams 
from  the  nanostructure.  This  allowed  us  to  make  the  first  observation  of  the  ballistic  heat  flow 
regime  at  nanoscale  interfaces  and  quantitatively  measure  the  strength  of  the  ballistic  effect. 
During  this  project  period,  we  have  focused  on  the  measurement  of  nanoscale  heat  transfer 


surrounding  even  smaller  heat  sources  (down  to  a  line  width  of  25  nm).  We  uncovered  a  new  and 
surprising  regime  of  nanoscale  thermal  transport  that  dominates  when  the  separation  between 
nanoscale  heat  sources  is  comparable  with  phonon  mean  free  paths  -  a  regime  that  has  not  been 
observed  or  predicted  to  date.  We  find  that,  very  counterintuitively,  nanoscale  heat  sources 
spaced  closely  together  can  cool  more  quickly  than  widely-spaced  heat  sources  of  the  same 
size!  This  most  recent  discovery  has  been  published  in  PNAS  -  Proceedings  of  the  National 
Academy  of  Sciences  of  the  United  States  of  America  (Vol.  112,  pp.4846-4851,  2015).  Such  a 
discovery  gives  unprecedented  method  to  characterize  the  differential  thermal  conductivity 
contributions  of  phonons  with  different  MFPs  and  to  benchmark  predictions  from  first-principles 
density  functional  theory  (DFT)  calculations. 

3).  We  have  made  significant  progresses  in  the  understanding  of  thermal  transport  in  emerging 
materials.  Working  with  our  collaborators  world-wide,  we  have  studied  thermal  transport  in 
many  nanostructured  materials,  including  hybrid  organic-inorganic  crystals,  thin  films  fabricated 
using  atomic  layer  deposition  (ALD)  and  molecular  layer  deposition  (MLD)  techniques,  and 
most  recently  the  emerging  two-dimensional  materials  including  graphene,  silicone,  and 
transition  metal  dichalcogenides. 


In  the  following,  we  highlight  some  of  the  afore-mentioned  achievements: 


Thrust  #1.  A  first-principles-based  atomistic  Green’s  function  study  on  phonon 
transmission  across  dissimilar  materials. 

In  recent  years,  the  atomistic  Green’s  function  (AGF)  approach  has  been  shown  to  be  an 
efficient  method  to  study  frequency-dependent  phonon  transport  across  interfaces  of  dissimilar 
materials.  The  phonon  transmission  is  highly  dependent  on  the  details  of  atomic  configuration 
and  interatomic  interaction  around  the  interface.  AGF  approach  has  been  applied  to  study  a  wide 
range  of  interfaces,  including  sharp  (smooth)  interfaces,  rough  interfaces,  interfaces  with 
vacancy  defects  and  alloyed  interfaces.  Often  the  empirical  potentials  are  used  in  AGF 
calculations  to  describe  the  interatomic  interactions.  The  interatomic  force  constants  from  the 
first-principles  calculations  have  also  been  integrated  with  AGF  method  for  studying  interfaces 
of  materials,  especially  when  the  empirical  potentials  are  not  readily  available. 

However,  it  would  be  quite  challenging  to  extract  the  interatomic  force  constants  of  realistic 
material  interfaces  from  the  first-principles  calculations.  Comparing  to  the  first-principles 
prediction  of  thermal  conductivity  of  bulk  crystals,  the  lattice  near  an  interface  of  two  dissimilar 
materials  is  likely  to  be  distorted  due  to  the  lattice  mismatch  and  the  difference  in  the  force  field 
experienced  by  the  atoms  in  the  interfacial  region.  A  large  supercell  would  be  required  for  the 
first-principles  calculations  to  capture  the  essential  characteristics  of  the  lattice-mismatched 
interface,  which  leads  to  severe  numerical  challenges.  The  mass  approximation  (MA)  that  was 
employed  to  calculate  the  thermal  conductivity  of  alloys  has  thus  been  used  to  extract 
interatomic  force  constants  of  the  interfacial  region  from  the  first-principles.  Under  the  MA,  the 
differences  in  the  lattice  constants  and  the  force  fields  between  the  two  dissimilar  materials  are 
ignored  and  only  the  difference  in  atomic  mass  is  taken  into  account.  It  was  recently  pointed  out 
that  the  MA  tends  to  overestimate  the  thermal  conductivity  due  to  the  neglect  of  the  local  force- 
field  difference.  Other  studies  showed  that  the  MA  under-predicts  the  contribution  of  high- 


frequency  phonons,  which  leads  to  a  lower  thermal  conductivity.  It  is  unclear  whether  MA  is  a 
reasonable  approximation  when  integrated  with  the  AGF  approach  to  calculate  the  phonon 
transmission  across  the  interfaces  of  dissimilar  materials  where  there  exist  the  differences  in  both 
lattice  constants  and  force  fields. 

In  this  work,  we  proposed  an  integrated  first-principles-based  AGF  approach  using  higher- 
order  force  constant  model  (HOFCM)  to  compute  phonon  transmission  across  interfaces  of 
dissimilar  materials,  with  the  details  shown  in  Figure  1. 


Interfacial  System 


Figure  1.  Numerical  simulation 
procedure  of  the  integrated 
DFT-AGF  approaches:  MA- 
AGF  and  HOFCM-AGF.  In  both 
approaches,  an  interfacial 
system  is  first  converted  to  an 
effective  virtual  crystal.  In  MA, 
the  harmonic  force  constants  are 
computed  from  DFT  calculations 
of  the  virtual  crystal,  and  then 
directly  used  as  the  inputs  for 
AGF  calculations.  In  HOFCM, 
the  interatomic  force  constants 
with  respect  to  both  the  atom 
displacements  and  atom  species 
are  extracted  from  DFT 
calculations.  Then  a  force-field 
using  these  interatomic  force 
constants  are  constructed,  which 
is  used  to  determine  the 
equilibrium  configuration  of  the 
interfacial  system  and  the 
corresponding  harmonic  force 
constants  that  are  used  in  AGF 
calculations. 


By  considering  the  lattice  distortion  and  the  species-dependent  local  force  field,  the  accuracy 
of  the  interatomic  force  constants  extracted  is  improved  compared  with  that  from  the  MA,  while 
the  required  computational  resources  are  significantly  less  severe  than  directly  extracting  force 
constants  from  the  standard  first-principles  calculations.  As  an  example,  we  conducted  the 
detailed  studies  of  frequency-dependent  phonon  transmission  across  Mg2Si/Mg2Sii-xSnx  interface, 
as  shown  in  Figure  2,  which  is  promising  thermoelectric  material  system  for  intermediate 
temperature  range  energy  conversion  applications.  We  systematically  study  how  the  lattice 
mismatch-induced  local  force  field  influences  phonon  transmission  across  the  interfaces  between 
dissimilar  materials,  which  are  inevitably  ignored  in  the  MA.  This  study  could  provide  guidance 
for  designing  nanostructured  materials  with  tunable  thermal  conductivity. 


(a) 


Left  reservoir 


Right  reservoir 
Mg2Si/Mg2Si1.xSnx  interface 


Figure  2.  (a)  The 

schematic  of  the 
simulation  system  of 
Mg2Si/Mg2Sii- 
xSnx/Mg2Si  structure, 
(b).  Phonon 

transmission  across  the 
Mg2Si/Mg2Sii- 
xSnx/Mg2Si  structure  as 
a  function  of  phonon 
frequency  for  different 
Sn  compositions  (x).  The 
length  of  the  alloy  layers 
is  fixed  at  10  unit  cells. 


Thrust  #2.  phonon  mean  free  paths  measurement  technique  using  extreme  ultraviolet 
laser-enabled  pump-and-probe  method. 

The  Fourier’s  law  of  heat  conduction  considers  heat  transport  as  a  diffusive  process  where 
energy  flow  is  driven  by  a  temperature  gradient.  However,  at  length  scales  smaller  than  the  mean 
free  path  for  the  energy  carriers,  heat  flow  becomes  ballistic  -  driven  by  direct  point-to-point 
transport  of  energy  quanta.  Although  past  experiments  have  seen  ballistic  effects  in  layered  thin 
films  or  nanowires,  non-Fourier  heat  dissipation  from  a  nanoscale  hot  spot  to  the  surrounding 
had  not  been  experimentally  observed  and  was  a  subject  of  some  uncertainty  although  it  has 
significant  relevance  in  heat  dissipation  of  electronic  transistors  In  an  earlier  work  published  in 
Nature  Materials  (Vol.  9,  pp.  26-30,  2010),  we  developed  a  pump-and-probe  characterization 
tool  using  table-top  femtosecond  extreme  ultraviolet  (EUV)  laser  to  study  nanoscale  thermal 
transport.  In  the  experiment,  an  array  of  nanowires  made  by  electron-beam  lithography  is  heated 
with  an  ultrafast  laser  pulse,  and  the  heat  flow  from  the  nanowire  array  into  a  bulk  substrate  is 
monitored  by  diffracting  ultrafast  coherent  EUV  beams  from  the  nanostructure.  This  allowed  us 
to  make  the  first  observation  of  the  ballistic  heat  flow  regime  at  nanoscale  interfaces  and 
quantitatively  measure  the  strength  of  the  ballistic  effect.  Such  an  experimental  discovery  of 
ballistic  phonon  transport  a  surrounding  nanoscale  heat  source  along  with  the  recent  advances  on 
computational  prediction  of  phonon  mean  free  path  distributions  using  the  first-principles-based 
density  functional  theory  have  inspired  a  few  research  groups  in  exploring  the  measurements  of 
phonon  mean  free  path  distributions  in  crystalline  materials  using  ultrafast  optical  wavelength 
laser-based  pump  and  probe  method  [Regner,  K.  T.  et  al.  Nat.  Commun.  4,  1640  (2013). 
Minnich,  A.  J.  et  al.  Phys.  Rev.  Lett.  107,  095901  (2011).  Johnson,  J.  A.  et  al.  Phys.  Rev.  Lett. 
110,025901  (2013).]. 

However,  the  basic  proof-of-principle  using  optical  wavelength  lasers  has  limited  to  the 
contribution  of  long-wavelength  (>  1  pm)  MFP  phonons  in  silicon.  During  this  project  period. 


we  have  focused  on  the  measurement  of  nanoscale  heat  transfer  surrounding  even  smaller  heat 
sources  (down  to  a  linewidth  of  25nm)  using  the  table-top  femtosecond  EUV  laser-enabled 
pump-and-probe  method.  We  uncovered  a  new  and  surprising  regime  of  nanoscale  thermal 
transport  that  dominates  when  the  separation  between  nanoscale  heat  sources  is  comparable  with 
phonon  mean  free  paths  -  a  regime  that  has  not  been  observed  or  predicted  to  date.  We  find  that, 
very  counterintuitively,  nanoscale  heat  sources  spaced  closely  together  can  cool  more  quickly 
than  widely-spaced  heat  sources  of  the  same  size!  This  most  recent  discovery  has  been  published 
in  PNAS  -  Proceedings  of  the  National  Academy  of  Sciences  of  the  United  States  of  America 
(Vol.  112,  pp.4846-4851,  2015).  Such  an  experimental  demonstration  of  the  regime  map  of 
nanoscale  heat  transport  gives  unprecedented  method  to  characterize  the  differential  thermal 
conductivity  contributions  of  phonons  with  different  MFPs  and  to  benchmark  predictions  from 
first-principles  density  functional  theory  (DFT)  calculations.  As  shown  in  Figure  3,  the 
experimental  data  of  the  size-dependent  effective  thermal  resistivity  (conductivity)  for  sapphire 
and  silicon  substrates  were  used  to  extract  the  phonon  mean  free  path  distributions  by  using  the 
suppression  function  which  identifies  the  contributions  of  phonon  modes  with  each  mean  free 
path  range.  To  extract  the  phonon  MFP  distributions,  the  experimental  data  were  fit  with  the 
weights  &(Ai)  assigned  to  multiple  bins  of  phonon  modes  which  give  their  average  relative 
contribution  to  the  differential  thermal  conductivity  (purple  shading).  Figure  3  shows  that  the 
experimental  data  on  cumulative  thermal  conductivity  and  phonon  mean  free  path  distributions 
agrees  well  with  DFT  calculation,  which  offers  the  first  experimental  benchmark  for  first- 
principle  calculations. 


Figure  3.  The  experimental 
data  on  effective  thermal 
resistivity  ( conductivity )  from 
nanoscale  heat  source  was 
used  to  extract  the 
cumulative  thermal 

conductivity  and  differential 
thermal  conductivity  due  to 
the  contributions  of  phonon 
with  different  mean  free 
paths.  The  phonon  mean  free 
path  distributions  measured 
agree  well  with  DFT 
calculation. 
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Thrust  #3.  Modeling  and  measurement  of  thermal  transport  in  emerging  materials. 

Working  with  our  collaborators  world-wide,  we  have  studied  thermal  transport  in  many 
nanostructured  materials,  including  hybrid  organic-inorganic  crystals,  thin  films  fabricated  using 
atomic  layer  deposition  (ALD)  and  molecular  layer  deposition  (MLD)  techniques,  and  most 
recently  the  emerging  two-dimensional  materials  including  graphene,  silicone,  and  transition 
metal  dichalcogenides. 

Two-dimensional  transition  metal  dichalcogenides  (TMDCs)  are  finding  promising 
electronic  and  optical  applications  due  to  their  unique  properties.  In  this  work,  we  systematically 
study  the  phonon  transport  and  thermal  conductivity  of  eight  semiconducting  single-layer 
TMDCs,  MX2  (M=Mo,  W,  Zr  and  Hf,  X=S  and  Se),  as  shown  in  Figure  4,  by  using  the  first- 
principles-driven  phonon  Boltzmann  transport  equation  approach.  The  validity  of  the  single¬ 
mode  relaxation  time  approximation  to  predict  the  thermal  conductivity  of  TMDCs  is  assessed 
by  comparing  the  results  with  the  iterative  solution  of  the  phonon  Boltzmann  transport  equation. 
A  very  high  thermal  conductivity  value  of  142  W/mK  was  found  in  single-layer  WS2.  The  large 
atomic  weight  difference  between  W  and  S  leads  to  a  very  large  phonon  bandgap  which  in  turn 
forbids  the  scattering  between  acoustic  and  optical  phonon  modes  and  thus  resulting  in  very  long 
phonon  relaxation  time,  as  shown  in  Figure  5. 


Figure  4.  The  thermal  conductivity  of  (a)  2H  and  (b)  IT  TMDC  monolayers  as  a  function  of 
temperature. 


Figure  5.  Phonon  dispersion  (a)  and  lifetime  (b)  of  MoS 2  and  WS2  calculated  from  the  first- 
principles  simulations. 
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Abstract 

The  Air  Force  is  increasingly  relying  on  the  benefits  of  novel  multifunctional  materials.  These  materials  are  being  pushed  to 
their  thermal  limits  in  progressively  more  demanding  roles.  Therefore,  there  is  a  strong  need  to  develop  modeling  and 
characterization  tools  for  better  understanding  of  thermal  transport  in  nanostructured  functional  materials.  The  objective  of  this 
project  is  to  develop  an  integrated  modeling  and  characterization  tool  to  study  phonon  transmission  and  generation  across 
engineered  and  strained  interfaces.  Such  a  tool  will  be  critical  in  understanding  the  fundamentals  of  energy  transport  and  in 
developing  structure-property  relations  of  functional  nanostructures  that  could  enable  novel  design  of  multi-functional 
nanostructures  with  superior  thermal  properties,  such  as  super-thermal  insulators,  high  efficiency  thermoelectric  materials, 
and  super-thermal  conductors. 


The  Innovation  and  Accomplishments  of  this  project  include: 


1) .  Developed  modeling  and  simulation  tools  for  thermal  and  thermoelectric  transport  processes  across  interfaces  and  in 
nanostructured  materials.  In  particular,  we  have  developed  a  novel  approach  to  predict  phonon  transport  properties  across 
dissimilar  interfaces.  Our  approach  is  based  on  non-equilibrium  Green's  function  theory  (NEGF)  which  provides  a 
wavelength-dependent  phonon  transmission  function.  This  approach  has  been  integrated  with  the  molecular  dynamics 
simulations  and  the  first-principles  simulations  to  predict  phonon  transmission  function  across  strained  and  engineered 
materials  interfaces.  The  obtained  wavelength-dependent  phonon  transmission  function  will  be  integrated  into  the  Boltzmann 
transport  model  to  study  and  design  large-scale  multidimensional  systems  with  embedded  nanostructures.  The  results  were 
published  in  Journal  of  Physics:  Condensed  Matter  (Vol.  24,  Art  #  1 55302,  2012),  Physical  Review  B  (Vol.  86,  Art  #054305, 
2012),  and  a  recently  accepted  paper  in  Physical  Review  B  (Xiaokun  Gu,  et  al,  Phonon  transmission  across  Mg2Si/Mg2Si1- 
xSnx  interface:  a  first-principles-based  atomistic  Green's  function  study,  http://arxiv.org/abs/1501 .04084,  May  2015) 

2) .  In  an  earlier  work  published  in  Nature  Materials  (Vol.  9,  pp.  26-30,  201 0),  we  developed  a  pump-and-probe 
characterization  tool  using  table-top  femtosecond  extreme  ultraviolet  (EUV)  laser  to  study  nanoscale  thermal  transport.  In  the 
experiment,  an  array  of  nanowires  made  by  electron-beam  lithography  is  heated  with  an  ultrafast  laser  pulse,  and  the  heat 
flow  from  the  nanowire  array  into  a  bulk  substrate  is  monitored  by  diffracting  ultrafast  coherent  EUV  beams  from  the 
nanostructure.  This  allowed  us  to  make  the  first  observation  of  the  ballistic  heat  flow  regime  at  nanoscale  interfaces  and 
quantitatively  measure  the  strength  of  the  ballistic  effect.  During  this  project  period,  we  have  focused  on  the  measurement  of 
nanoscale  heat  transfer  surrounding  even  smaller  heat  sources  (down  to  a  linewidth  of  25nm).  We  uncovered  a  new  and 
surprising  regime  of  nanoscale  thermal  transport  that  dominates  when  the  separation  between  nanoscale  heat  sources  is 
comparable  with  phonon  mean  free  paths  -  a  regime  that  has  not  been  observed  or  predicted  to  date.  We  find  that,  very 
counterintuitively,  nanoscale  heat  sources  spaced  closely  together  can  cool  more  quickly  than  widely-spaced  heat  sources  of 
the  same  size!  This  most  recent  discovery  has  been  published  in  PNAS  -  Proceedings  of  the  National  Academy  of  Sciences  of 
the  United  States  of  America  (Vol.  1 1 2,  pp.4846-4851 , 2015).  Such  a  discovery  gives  unprecedented  method  to  characterize 
the  differential  thermal  conductivity  contributions  of  phonons  with  different  MFPs  and  to  benchmark  predictions  from  first- 
principles  density  functional  theory  (DFT)  calculations. 

3) .  We  have  made  significant  progresses  in  the  understanding  of  thermal  transport  in  emerging  materials.  Working  with  our 
collaborators  world-wide,  we  have  studied  thermal  transport  in  many  nanostructured  materials,  including  hybrid  organic- 
inorganic  crystals,  thin  films  fabricated  using  atomic  layer  deposition  (ALD)  and  molecular  layer  deposition  (MLD)  techniques, 
and  most  recently  the  emerging  two-dimensional  materials  including  graphene,  silicone,  and  transition  metal 
dichalcogenides. 
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